LTK mutations responsible for resistance to lorlatinib in non-small cell lung cancer harboring CLIP1-LTK fusion.

The CLIP1-LTK fusion was recently discovered as a novel oncogenic driver in non-small cell lung cancer (NSCLC). Lorlatinib, a third-generation ALK inhibitor, exhibited a dramatic clinical response in a NSCLC patient harboring CLIP1-LTK fusion. However, it is expected that acquired resistance will inevitably develop, particularly by LTK mutations, as observed in NSCLC induced by oncogenic tyrosine kinases treated with corresponding tyrosine kinase inhibitors (TKIs). In this study, we evaluate eight LTK mutations corresponding to ALK mutations that lead to on-target resistance to lorlatinib. All LTK mutations show resistance to lorlatinib with the L650F mutation being the highest. In vitro and in vivo analyses demonstrate that gilteritinib can overcome the L650F-mediated resistance to lorlatinib. In silico analysis suggests that introduction of the L650F mutation may attenuate lorlatinib-LTK binding. Our study provides preclinical evaluations of potential on-target resistance mutations to lorlatinib, and a novel strategy to overcome the resistance.

Integrated approach to functional analysis of an ERBB2 variant of unknown significance detected by a cancer gene panel test.

PURPOSE: Dealing with variants of unknown significance (VUS) is an important issue in the clinical application of NGS-based cancer gene panel tests. We detected a novel ERBB2 extracellular domain VUS, c.1157A > G p.(E401G), in a cancer gene panel test. Since the mechanisms of activation by ERBB2 extracellular domain (ECD) variants are not fully understood, we aimed to clarify those mechanisms and the biological functions of ERBB2 E401G. METHODS: ERBB2 E401G was selected as VUS for analysis because multiple software tools predicted its pathogenicity. We prepared ERBB2 expression vectors with the E401G variant as well as vectors with S310F and E321G, which are known to be activating mutations. On the basis of wild-type ERBB2 or mutant ERBB2 expression in cell lines without ERBB2 amplification or variants, we evaluated the phosphorylation of human epidermal growth factor receptor 2 and related proteins, and investigated with molecular dynamics (MD) simulation the mechanisms conferred by the variants. The biological effects of ERBB2 E401G were also investigated, both in vitro and in vivo. RESULTS: We found that ERBB2 E401G enhances C-terminal phosphorylation in a way similar to S310F. MD simulation analysis revealed that these variants maintain the stability of the EGFR-HER2 heterodimer in a ligand-independent manner. Moreover, ERBB2 E401G-transduced cells showed an increased invasive capacity in vitro and an increased tumor growth capacity in vivo. CONCLUSION: Our results provide important information on the activating mechanisms of ERBB2 extracellular domain (ECD) variants and illustrate a model workflow integrating wet and dry bench processes for the analysis of VUS detected with cancer gene panel tests.

Structure-Based de Novo Molecular Generator Combined with Artificial Intelligence and Docking Simulations.

Recently, molecular generation models based on deep learning have attracted significant attention in drug discovery. However, most existing molecular generation models have serious limitations in the context of drug design wherein they do not sufficiently consider the effect of the three-dimensional (3D) structure of the target protein in the generation process. In this study, we developed a new deep learning-based molecular generator, SBMolGen, that integrates a recurrent neural network, a Monte Carlo tree search, and docking simulations. The results of an evaluation using four target proteins (two kinases and two G protein-coupled receptors) showed that the generated molecules had a better binding affinity score (docking score) than the known active compounds, and the generated molecules possessed a broader chemical space distribution. SBMolGen not only generates novel binding active molecules but also presents 3D docking poses with target proteins, which will be useful in subsequent drug design. The code is available at https://github.com/clinfo/SBMolGen.

Exploring ligand binding pathways on proteins using hypersound-accelerated molecular dynamics.

Capturing the dynamic processes of biomolecular systems in atomistic detail remains difficult despite recent experimental advances. Although molecular dynamics (MD) techniques enable atomic-level observations, simulations of "slow" biomolecular processes (with timescales longer than submilliseconds) are challenging because of current computer speed limitations. Therefore, we developed a method to accelerate MD simulations by high-frequency ultrasound perturbation. The binding events between the protein CDK2 and its small-molecule inhibitors were nearly undetectable in 100-ns conventional MD, but the method successfully accelerated their slow binding rates by up to 10-20 times. Hypersound-accelerated MD simulations revealed a variety of microscopic kinetic features of the inhibitors on the protein surface, such as the existence of different binding pathways to the active site. Moreover, the simulations allowed the estimation of the corresponding kinetic parameters and exploring other druggable pockets. This method can thus provide deeper insight into the microscopic interactions controlling biomolecular processes.

E487K-Induced Disorder in Functionally Relevant Dynamics of Mitochondrial Aldehyde Dehydrogenase 2.

Mitochondrial aldehyde dehydrogenase 2 (ALDH2), which is a homotetramer assembled by two equivalent dimers, is an important enzyme that metabolizes ethanol-derived acetaldehyde to acetate in a coenzyme-dependent manner. The highly reactive acetaldehyde exhibits a toxic effect, indicating that the proper functioning of ALDH2 is essential to counteract aldehyde-associated diseases. It is known that the catalytic activity of ALDH2 is drastically impaired by a frequently observed mutation, E487K, in a dominant fashion. However, the molecular basis of the inactivation mechanism is elusive because of the complex nature of the dynamic behavior. Here, we performed microsecond-timescale molecular dynamics simulations of the proteins complexed with coenzymes. The E487K mutation elevated the conformational heterogeneity of the dimer interfaces, which are relatively distal from the substituted residue. Dynamic network analyses showed that Glu487 and the dimer interface were dynamically communicated, and the dynamic community further spanned throughout all of the subunits in the wild-type; however, this network was completely rearranged by the E487K mutation. The perturbation of the dynamic properties led to alterations of the global conformational motions and destabilization of the coenzyme binding required for receiving a proton from the catalytic nucleophile. The insights into the dynamic behavior of the dominant negative mutant in this work will provide clues to restore its function.

Landscape and function of multiple mutations within individual oncogenes.

Sporadic reports have described cancer cases in which multiple driver mutations (MMs) occur in the same oncogene1,2. However, the overall landscape and relevance of MMs remain elusive. Here we carried out a pan-cancer analysis of 60,954 cancer samples, and identified 14 pan-cancer and 6 cancer-type-specific oncogenes in which MMs occur more frequently than expected: 9% of samples with at least one mutation in these genes harboured MMs. In various oncogenes, MMs are preferentially present in cis and show markedly different mutational patterns compared with single mutations in terms of type (missense mutations versus in-frame indels), position and amino-acid substitution, suggesting a cis-acting effect on mutational selection. MMs show an overrepresentation of functionally weak, infrequent mutations, which confer enhanced oncogenicity in combination. Cells with MMs in the PIK3CA and NOTCH1 genes exhibit stronger dependencies on the mutated genes themselves, enhanced downstream signalling activation and/or greater sensitivity to inhibitory drugs than those with single mutations. Together oncogenic MMs are a relatively common driver event, providing the underlying mechanism for clonal selection of suboptimal mutations that are individually rare but collectively account for a substantial proportion of oncogenic mutations.

Coarse-Grained Diffraction Template Matching Model to Retrieve Multiconformational Models for Biomolecule Structures from Noisy Diffraction Patterns.

Biomolecular imaging using X-ray free-electron lasers (XFELs) has been successfully applied to serial femtosecond crystallography. However, the application of single-particle analysis for structure determination using XFELs with 100 nm or smaller biomolecules has two practical problems: the incomplete diffraction data sets for reconstructing 3D assembled structures and the heterogeneous conformational states of samples. A new diffraction template matching method is thus presented here to retrieve a plausible 3D structural model based on single noisy target diffraction patterns, assuming candidate structures. Two concepts are introduced here: prompt candidate diffraction, generated by enhanced sampled coarse-grain (CG) candidate structures, and efficient molecular orientation searching for matching based on Bayesian optimization. A CG model-based diffraction-matching protocol is proposed that achieves a 100-fold speed increase compared to exhaustive diffraction matching using an all-atom model. The conditions that enable multiconformational analysis were also investigated by simulated diffraction data for various conformational states of chromatin and ribosomes. The proposed method can enable multiconformational analysis, with a structural resolution of at least 20 Å for 270-800 Å flexible biomolecules, in experimental single-particle structure analyses that employ XFELs.

Improvement in predicting drug sensitivity changes associated with protein mutations using a molecular dynamics based alchemical mutation method.

While molecular-targeted drugs have demonstrated strong therapeutic efficacy against diverse diseases such as cancer and infection, the appearance of drug resistance associated with genetic variations in individual patients or pathogens has severely limited their clinical efficacy. Therefore, precision medicine approaches based on the personal genomic background provide promising strategies to enhance the effectiveness of molecular-targeted therapies. However, identifying drug resistance mutations in individuals by combining DNA sequencing and in vitro analyses is generally time consuming and costly. In contrast, in silico computation of protein-drug binding free energies allows for the rapid prediction of drug sensitivity changes associated with specific genetic mutations. Although conventional alchemical free energy computation methods have been used to quantify mutation-induced drug sensitivity changes in some protein targets, these methods are often adversely affected by free energy convergence. In this paper, we demonstrate significant improvements in prediction performance and free energy convergence by employing an alchemical mutation protocol, MutationFEP, which directly estimates binding free energy differences associated with protein mutations in three types of a protein and drug system. The superior performance of MutationFEP appears to be attributable to its more-moderate perturbation scheme. Therefore, this study provides a deeper level of insight into computer-assisted precision medicine.

Rotation Mechanism of Molecular Motor V1-ATPase Studied by Multiscale Molecular Dynamics Simulation.

Enterococcus hirae V1-ATPase is a molecular motor composed of the A3B3 hexamer ring and the central stalk. In association with ATP hydrolysis, three catalytic AB pairs in the A3B3 ring undergo conformational changes, which lead to a 120° rotation of the central stalk. To understand how the conformational changes of three catalytic pairs induce the 120° rotation of the central stalk, we performed multiscale molecular dynamics (MD) simulations in which coarse-grained and all-atom MD simulations were combined using a fluctuation matching methodology. During the rotation, a catalytic AB pair spontaneously adopted an intermediate conformation, which was not included in the initial inputs of the simulations and was essentially close to the "bindable-like" structure observed in a recently solved crystal structure. Furthermore, the creation of a space between the bindable-like and tight pairs was required for the central stalk to rotate without steric hindrance. These cooperative rearrangements of the three catalytic pairs are crucial for the rotation of the central stalk.